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We present some consequences of non-anomalous propagation requirements on vari- 
ous massless fields. Among the models of nonlinear electrodynamics we show that only 



> 

(T) ■ Maxwell and Born-Infeld also obey duality invariance. Separately we show that, for 

in 

\ actions depending only on the F^ v invariant, the permitted models have L ~ vl + F 2 - 

o 

We also characterize acceptable vector-scalar systems. Finally we find that wide classes 
of gravity models share with Einstein the null nature of their characteristic surfaces. 

<D 

' The nonlinear electrodynamics of Born-Infeld (BI) has been the subject of frequent re- 

> : 

vivals, not least because it enjoys two quite separate properties, shared with Maxwell theory. The 
first is that its excitations propagate without the shocks common to generic nonlinear models 
||, |[] j the second is duality invariance ||. Here, we want to complete this subject in several 
ways. Our primary result will be that imposition of both duality and good propagation singles 
out BI and Maxwell, without even requiring the solutions to reduce to Maxwell in the weak field 
limit (something that is used to select BI in the derivations of each separate demand). We will 
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then show separately that for actions depending only on the Maxwell invariant a = \Ff LV (rather 
than on both a and the other invariant, (3 = \e^ uaT F^ v F aT ^ allowed in D=4) correct propagation 
implies the form L ~ y/1 + a, which shares a square root (but not a determinant) form with BI. 
In this connection, we will also discuss the scalar analogs of BI and characterize systems involving 
both scalars and vectors. Criteria for physical propagation of hyperbolic systems can be derived in 
many ways; in a separate work || we will consider the "complete exceptional" (CE) approach |?J 
to these questions, in more detail. Since this method was useful in some of our derivations below, 
we will use "CE" to mean propagation without emerging shocks, with characteristic flow remaining 
parallel along the waves (but waves need not travel with the same speed in all directions.) 

In the following we will be concerned with systems of partial differential equations which arise 
as the variational equations of a relativistic Lagrangian theory. Thus they will be quasilinear (linear 
in highest derivatives) and the coefficient functions will not explicitly depend on the coordinates. 
Without loss of generality, they can be reduced to a set of equations first order in derivatives. 
For definiteness then, with u an iV-vector of fields, A an NxN matrix and B an iV-vector (both 
arbitrary (smooth enough) functions of u), the equations of interest are 

^(«)a M «+B(«) = 0. (i) 

Of course, the theory of such equations in arbitrary dimensions is quite difficult, but we will be 
interested in the evolution of the spatial boundary of a wave propagating into some given vacuum 
state. So, with u some smooth (say at least C 1 ) solution, at some initial time we have some spatial 
region outside of which the "state" is the "vacuum solution" u, and across the boundary surface 
the full solution u is continuous but its first derivative may not be. We now consider the evolution 
of this initial "wavefront". 

One discussion of this situation has been developed (|| and references therein) which we will 
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follow here. Let the hypersurface S, specified by 

p(x) = 0, (2) 

denote the surface of evolution of the initial wavefront; i.e., the initial wavefront is the spatial 
surface (p(x, 0) = 0. Since the field u is continuous, only the normal derivative will be allowed to 
be discontinuous. Choosing a local coordinate system x^ = {ip,ip l \ we can define the "first order 
discontinuity" in a given quantity / to be 

df 



6if = 
where 



dip 



(3) 



[X] = X\^ =0+ - X|^ =0 - (4) 

and we will sometimes write [X] = SoX. 

Taking the discontinuity of (|l|) we obtain (ip^ = d^tp) 

(A"ip t ,)S 1 u = 0, (5) 

where the matrix A = A(u). Hence, since 5\u ^ 0, we see that S must be a characteristic surface; 
i.e., the characteristic equation 

H{x, <p„) = det (Afyp) = (6) 

must hold on S, where H is homogeneous of order N in ip^. The characteristic curves, which solve 

dx* _ 8H 
ds dp^ ' 

are clearly tangential to the characteristic surface. 

From (||), it further follows that 5iu may be expanded in a basis of Ker (A^p^,), the coefficients 
of this expansion being called the "coefficients of discontinuity". In general, the coefficients of 



discontinuity evolve according to a nonlinear differential equation. Following the discussion in 
[]|], based on ||, it can be shown that when we ensure that the characteristic curves do not 
intersect locally (and thus "shock" singularities do not develop), the evolution of the equations 
of discontinuity is correspondingly linear. We will impose this condition as our specification of 
"nonanomalous evolution". It is the CE condition mentioned earlier, and can be imposed as the 
condition that, on the characteristic surface H = 0, we have Q 

5 H = 0. (7) 

We now apply these ideas to the systems of interest, without however showing any details 
of the underlying CE derivation. For general nonlinear gauge invariant actions in D = 4 that 
depend only on field strengths but not explicit derivatives, the Lagrangian L(a,f3) must obey the 
two equations (subscripts indicate partial differentiation) 

-L a (4 L aa - Lpp) + 2 a [L aa Lpp - (L a/3 ) 2 ] = , (7a) 

—L a Lap + (3 [L aa Lpp — (L a/3 ) 2 ] = . (76) 

in order for its excitations to be CE and for light to travel according to only one dispersion law, 
i.e. no birefringence ||, Q. To obtain them, however, it was necessary to assume nontrivial j5 
dependence; see below. Quite independently the duality invariance requirement is given by H 

{L a f-^L a L p -\{L p f= l - . (8) 

The simultaneous solution of (7) and (|8|) can be obtained as follows. 

We first change to variables (a, 6) with a = a + b, 1 = — ab (these new variables have the 
fundamental advantage that they factorize L^i)', then (7a), (76) and (||) transform into 

-4a6(a + 6)(a - bf[L aa L bb - (L ab ) 2 } + 4a6(a 2 - b 2 )L ab {L a - L b ) + 8a6(6 - a)L ab {bL b - aL a ) 
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-2(a - b) 2 (a + b)[aL b L aa + bL a L bb ] + (4a6 + (a + bf){L a - L b ){aL a - bL b ) = 0, (9a) 
2ab(a - b) 2 [L aa L bb - {L ab ) 2 } + (a + b)(L a - L b )(bL b - aL a ) + (a - b) 2 [aL b L aa + bL a L bb ] 

-{a 2 - b 2 )L ab (aL a - bL b ) + 2ab(b - a)L ab {L a - L b ) = 0, (96) 

L a L b = \ , (10) 

respectively. Note the symmetry in (a, 6) of each equation. Multiplying (96) by 2(a + 6) and adding 
the result to (9a), one finds 

(a - 6) 2 (aL a - 6L fe )[2(6 - a)L ab - (L a - L b )] = . (11) 



Vanishing of the (aL a — bL b ) factor is not a useful solution of ( [TT] ) because it would imply L = 
L{ab) = L(f3 2 ), for which (||) implies L ~ (3, a total divergence. Hence we must impose 

2(6 - a)L ab - {L a - L b ) = . (12) 

Substituting this back into (9a) gives the simplified form 

2ab[L aa L bb - (L ab ) 2 } + [a{L a L b ) a + b(L a L b ) b ] = . (13) 



So now instead of (9) and pQ), we can study (fjjh, (lq) and (^0[). But actually it is easy to see that 



( |l3| ) follows from (10), since (10) is a first integral of the Monge- Ampere factor [L aa L bb — (L, 



ab) 2 ] 



i.e. this factor vanishes as a consequence of (|10D. Thus we are left with the system (|l0D, ( |12| 



Substituting for L b in (12), it can be written as (a — 6) 2 { ^^^ ) = 1, with first integral 



L a = -^/l + (a-6)/(a) . (14) 



4 An alternative derivation inserts the general parametric solution H of ( |lo[ ) into (p"2|). 
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But by ©, 4(L a ) 2 = L a /L b = 1 + (a - 6)/(a) which gives L b = \[l + (a - ^/(a)^ 1 / 2 , whose 
integral is 



+ + . (15a) 

Given the (a, 6) symmetry of the equations, the above procedure based on L b instead of L a gives 
the corresponding form 

L = --— y/l + {b- a)k(b) + m(b) . (156) 



Consistency with (10) demands that -^/l + (a — 6)/(a)^/l + (b — a)k(b) = 1, which implies k(b) 
i + (/ib)/( a ) • Substituting for /c(6) in (156) then shows that 



L = -j^-^l + (a-b)f(a)+m(b) (16) 



which implies h(a) = m(b) = d = const. Differentiating (|T(j) with respect to a and comparing the 
result with (|l4j) one finally finds that (2 + (a — b)f(a))(f — f 2 ) = 0, which when integrated gives (for 
s an integration constant) f(a) = (or / = trivially, the Maxwell case.) Finally then, renaming 
constants, L = + ac)(l + be) + d (and L = — |(a + b) for / = 0). Rewriting these using 

(a,/3), we find using allowed rescalings that L = —\ct, Maxwell, and L = 1 — y/l + a — f3 2 = 1 — 



det^j, + Ffj, u ], Born-Infeld, are the only possible solutions. This demonstrates the uniqueness 
of BI and Maxwell as the simultaneously CE and duality invariant electrodynamics. The CE 
requirement (7) alone permits additional solutions, such as L = a/ f3, that are not duality invariant. 
[Conceivably, requiring power series expandability in (a, (3) might restrict the solutions of (7) to 
Maxwell and BI.] 

As was mentioned earlier, the CE criteria (7) are only correct for nontrivial /5-dependence; for 
pure L(a) they just imply V = const., namely Maxwell. Instead, if one goes back to the complete 
CE requirements, they imply, for pure L(a) 

L'L'" - 3{L") 2 = . (17) 



The solution, apart from Maxwell (L" = 0) is 

L(a) = k + (d + ca)^ . (18) 

[We remark that in D = 3, where a is the only invariant, this is also the CE result, there also 
Vl + oi = \J — det[r)^ u + F^ v ] with the BI determinant form. In D = 2 there is of course no 
propagation for any L(a) and correspondingly no restrictions are imposed.] The above BI form 
is very analogous to that obtained for a scalar field 0|. There, (in any dimension) for L(z), 
z = ^(d ll 4>)(d^(f>) being the only invariant (in first derivatives), we find the same requirement (|l7|). 
Hence we obtain, apart from the free scalar solution L = — \z, the same square root solution. 
Amusingly, this one can be put into Bl-like form, since (rescaling </>) 

VT+2z~ = + = ^-detfo^ + ^l (19) 

where (j)n = d^(j) denotes the field strength. If one combines Maxwell and the (neutral) scalar into an 
action L(a,P,z), then the CE conditions further require L za = = L z p, reducing the Lagrangian 
to the noninteracting L(q, (3)+L(z) form. Having the "fully" BI form J ' — detfry^ + F^ v + (j)^^ i n 
mind, one can actually show more generally that there are no CE actions with nontrivial dependence 
on the other possible variable y = ^{F jiu (f) l/ ) 2 . 

Finally, we turn to gravitation. For Einstein's gravity in vacuum, as well as the linearized 
theory, the gravitational waves are CE, the characteristic surfaces describing discontinuities being 
null (see e.g. |l(|). It can be shown that this result holds for any D > 4. [For D = 3, there is of 
course no propagation and no restrictions are imposed.] One can further look at pure gravitational 
actions of the form pR^ — qR 2 in D = 4 and f(R) in D > 3 and show that the same conclusion 
remains unchanged. 

To reduce these theories to a first order system would be inconvenient, but is fortunately 
made unnecessary by a simple extension of the previous discussion. Clearly, if we rebuilt the 
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original higher order equations from the set (||), we would simply have the situation that all the 
derivatives of the field are assumed continuous except the highest one. Thus for quasilinear systems 
higher order, say q, in derivatives, we define 



S r f 



d r f 



d<p r 



and will consider the case that 



S g u / ; S r u = , < r < q . 



Notice that 



(5 r 3 M = d^f 5 r+ i 



Let us first sketch the Einstein case to establish notation. Considering a second order dis- 
continuity in the metric across some characteristic surface <p = 0, 829^ = we have (ip^ = d^ip) 



and 



S R = grtfoRp,) = Wnp, ~ <P^k u u 



which implies for 



5nG 



SqG^u — So(R^ u — -g^R) — SqR^ u — -g^SoR — 



^M^A7T A v + VWA7T A fj, ~ V^vK X X ~ V ? A^ A ^ ~ V^V^ctt ~ ^a^^ T ) 



. 

(20) 
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In the harmonic gauge g^T a ^ = 0, one finds that its first discontinuity implies 

2tt^^ - ^ ^ v = (21) 
Multiplying this by g va tp T + 5W<Ax, one gets 

^^avt A u + tPuP\-K X n - (p^uTr x x = , (22) 
whereas contracting by ip u , one finds 

^%, = ^W^7T% . (23) 



Using (E2J) and (23) in (p0[), one ends up with 



SqG^ = ^xf X ^ u + ^¥>A¥>V a) = 



Hence taking the trace 

<5oG" „ = ^±^v\^ a = 

The discontinuity in g^ u is arbitrary, hence ir u a ^ 0, which implies that <^a¥ :,A = 0. This tells 

that the characteristic surfaces are null: the discontinuities travel with the speed of light in all 

directions. The same holds for the linearized version of the theory as well of course. 

For generic quadratic Lagrangians (pRui/R^ —qR 2 )^/—g in D = 4, using similar steps (writing 

the field equations, choosing harmonic gauge as before and utilizing the identities (^), (|23|)) one 

finds that (Q = ip x tfx , ^ = tt A a) 

/ 1 pip \ 

Q[^(P- 2q)ip^ip u 7r - -Qtt^ v - ^9^(^ ~ 2q)Qir 1=0. (24) 

Taking the trace, one gets Q 2 ir(p — 3q) = 0. (The choice p = 3q corresponds to Weyl-tensor 
squared; the scalar degree of freedom is absent.) For p = 3q, (24) becomes 



1 3 1 

qQi-^PvVvK - -Qn^u + -g^Qir) = 



Since tt^ is arbitrary, we see that again Q = 0, as in Einstein, so Q = characterizes both Einstein 
and the quadratic action. 

For actions f{R) y / zr g in D > 4, the field equations are 

= R^f - \g^f + (^V CT V ff - V„V u )f = . 

Hence the order of highest derivatives is four. Following similar steps by taking b^g^v = t^uvi we 
find the same expressions for 5^T X ^ v and 52R 11U as for 5\T X ^ and SqR^u in the Einstein case. 
Using these, we get 

SaEfa, = (Qg^u - ^P^u)(^ a ^ T ^aT - Qn)f" = . 
Finally, in the harmonic gauge with identity ( ^3|) and taking the trace, one gets 

Here too Q = is the only solution, and so for a wide class of gravitational actions the propagation 
obeys the Einstein behavior. 

This work was supported in part by NSF, under grant no PHY-9315811. 
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